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ABSTRACT 


This report summarizes the progress under NASA Grant NSG 1213 
from July 1, 1975 to November 30, 1977. Section I provides a pre- 
sentation of the problem under consideration. In Section II, we 
consider the effects of system parameter variations on the over-all 
system stability. Section III contains results from the application 
of numerical minimization to reduce system coupling. In Section V, 
some future research directions are outlined. 
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I, Problei description 

Consider a i;iven linear time-invarient system described by the 
following equations; 


x(t) - Ax(t) + Bu(t) + LWg(t) 

(la) 

y„(t) - Cx(t) + Wg(t) 

(lb) 

y^j(t) - Dx(t) 

(Ic) 


where x(t) is the state vector, u(t) is the input vector, W|r(t) iS 

♦ 

the gust disturbance, ym(t) is the measurable output, y^Ci) is the 
output vector to be controlled. The gust Wg(t) is assumed to 
satisfy the following equation: 

wg(t) - Fj^Wg(t) + F2‘n(t). 

T) and ws are white noises. 

In the above system, we assume that the matrices A B, L, F]^ 
and F2 may vary under different operating conditions. Let Aq and 
Bq denote the nominal values of A and B at soiao chosen operating 
condition. 

Since the full state vector x(t) is not directly accessible, 
the following full-order state observer is used to estimate x(t); 

x(t) « (AQ-GC)x(t) + BQu(t) + GLCx(t)+Wg(t)] (2) 

where G matrix is chosen to assign a set of desired poles for (Aq-GC) 
in order to guarantee x(t)-»x(t) as t-» ». 

The above estimated state x(t) is used in the follcAving feedback 
control law; 


u(t) - ,Fv(t) +Kx(t) 
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vhere v(t) is the external command input. B/ usin^; decoupling 
theory, one can find appropriate matrices F and K such that the 
transfer function from v to is diagonal. 

The complete system is described by the lollowing equations: 


x(t) « Ax(t) + B ws 

T\ 


(3a) 


y^(t) - Cx(t) 


(3b) 


‘a BK L 

GC Aq-GC+BqK 0 

0 0 Fi 


BqF G 0 
0 0 F2 


C «[ D 0 0 ] 


The following diagram shows the complete structure of the whole 


Pla nt 

A - Ax + Bu + Lw 
ym “ Cx + Ws 
Yc " Dx 


B—i/ 

Tjd 


OBservet 

x-’(Ao*-GC)x+HqU 
+ G(Cx+Ws) 


u 
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The aircraft model studied is based on the report by Miller, 

Deal and ChampineLO. The system input consists of throttle deflection, 
horizontal~tail deflection and flap deflection. The system state 
consists of pitch anp:le, pitch rate, angle of attack and forward 
speed. The system outputs to be controlled are flight path angle, 
pitch angle and forward speed. The measurable outputs are pitch rate 
and forward speed. For different values of angle of atcack, the 
system parameters A and B vary. We choose a = 10^ as the trim flight 
condition and its parameters as the nominal A and B . From the 

o o 

decoupling theory, the decoupled system transfer function has the 
following form: 




^1 0 0 
s - on 


Vy 

3 

= 

0 ^2 0 

SS-021S-022 


<> 

CD 

r- — 
c> 
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0 0 ^3 . 

S-0 31 




where Xf , X2 , \3 , , 031 - ^22 * *^3l arbitrarily assigned 

by choosing appropriate feedback and feedforward matrices F and K. 

The external disturbance » on the angle of attack tr is due to 
the vertical wind of gust w , 



u is the forward velocity and Wg is assumed to have Dryden power 
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spectrum density of the form: 


«(n) 


„ iw 1+3 (twn)® 

V n TTTlZwnPT^ 


Iw is the scale of turbulence, is the variance of ihvi j?ust 
velocity. Hence be modeled as the output of a secon<l-or'der 

filter 


T(s) 



1 + /3 -tws 
( 1 + tws ) ® 


driven by white noise T). The sensor noises iu the measurement of 

pitch rate and forward speed are assumed to bo white noises w , and 

si 

Wg 2 respectively. 

The complete system has the following configuration : 
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“■ £gra"eter variations and s ystems 

The design of observer end feedbeck control lew Is besed on the nomine! 
system peremeters end corresponding to 10 degrees of engle of etteck. 
we heve essumed thet the engle of etteck veries between 5 degrees end 15 
Agrees. Different sets of poles heve been assigned to the observer end the 
system besed on the nominal system parameters end B . It turns out 
that the overall system poles ere very sensitive to system" parameter variations 
in Figure 1 and Figure 2. the overall system root locus ere plotted as the angle 
of attack varies from 5 degrees to 16 degrees. For the case in which the nominal 
system poles are -8. -5. -5, and -10. the system remains stable in the whole 
range of variations of « . However, if one assigns -8. -9 a 7j end -10 as the 
nominal system poles, the overall system becomes unstable as a approaches 5 
degrees. In the former case, the time history of the system response are shown 
in Figure 3-6. It can be seen that system parameter variations have caused 
Significant coupling among system inputs end outputs. Hence the straightforward 

application of observer theory and decoupling theory to flight control system 
nay not always be practical. 
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III. Reduction of coupling by using numerical minimization algorithm 
In this section, we apply Box Complex Algorithm [21 to the 
adjustment of feedback and feedforward control law (F, K) in order 
to reduce the couplinj: of the system under painmeCer variations. 

Two different cases were tried; each with a diL'f’orent function to 
be minimized and different explicit and implicit constraints. 


Case 1 


The function to be minimized is 

3 



i<*5°,10° 15' 


L 

, j”l 
3 


10 

> 

w=0 




2 


which is the sum of the magnitude squared of tl;e off-diagonal 
elements of the over all system transfer function when alpha is 
5°, 10° and 15°. 

There are three constraints imposed when searching for optimal 

(F, K), 

1) All elements of F and K, (totalling 21 in our ctfie) , can 
differ from the initial value of F and K (computed by 
Gilbert's decoupling program) by no more than ± 15.0. This 
provides an explicit feasible region lor the search of F and 
K. 

2) The determinat of F must be greater than certaiii value. 

This prevents the trival solution where I’ becomes a zero 


matrix. 
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3) The new pair (F.K) must produce a stable over-all system. 

This requirement is {piaranteed by using Kouth criterion. 
The results are summarized as follows: 
l3efore optimization - Fqj^ “ 13.52. 

After optimization - F^^ * 0.94. 


Let 


3 

10 


'T 



F 

'dg 

- L 

i 

u 

L 

»(=5°,10°,15° 


w=0 


be a measure of the magnitude of the diagonal olo:nents of tlie over 
all transfer function. Then we have the following: 

Before optimization - Fjj^ » 57.40. 

After optimization - Fj^^ - 33.90. 

In the following table, we give details on i;he mag.iitudes of 
diagonal and off-diagonal elements of the over-all system transfer 
function when * 5°, 10° and 15°. 



ALPHA 

(degrees) 

Diagonal 

Elements 

Of f-Di agonal 
Elements 

I" 

Before 

Opt. 

5 

21.54 

9.07 

10 

18.86 

5.5 X lO"® 


15 

17.01 

4.45 

After 
Opt . 

5 

10.59 

0.787 

10 

11.76 

0.0213 


15 

11.54 

0. 130 
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Althou{;l\ there seems to be a substantial improvement in the » 
reduction of coiiplinp: nt'ter optimization procedures, however, time i 
histories of the optimized system show that Ihore is not a Kif'nili- 
cant dil’ferencc. In fact, some steady state values of the optimized 
system actually increased as shown below; 

Let H(s) be the transfer function of the decoupled system, i.e.. 



Each column has t>een normalized so that the diagonal elements 
are 100. This result suggests that minimization of the steady state 
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vali|«s of the coupling iiight improve the system performance. 
This has been carried out as follows: 


Case 2 


The function to be minimized is 

3 


SOD 


L L 

U -5°, 10°, 15° 


h ,( 0 )| 
i.l I 


which is the sum of the square of the steady-state value of the off- 
diagonal elements of transfer function, for alpha equals 5° , 10° and 
15°. 


The following constraints are imposed for the optimization 
algorithm: 

1) All elements of F and K can differ from the initial 
values of F and K (computed by Gilbert's decoupling 
program for the nominal system, i.e, when e( » 10 °) by no 
more than * 

2) The diagonal elements of H(0) , i.e, h^^(O), 1122 ( 0 ) and 
1123 ( 0 ) are constrained to be greater than certain positive 
values. 

3) Every pair (F,K) should result in a stable system. 

The results are as follows: 

Before optimization - - 7.85 X lo”^ 

After optimization - - 6,73 X lo”^. 


F 

'SDG 


I 

L. 

»a'-5°,10° ,15° 


I 


fc*«l 



(0) 


Let 


2 



16 


Before optimization - - 8.95 

After optimization - ^gjjQ * 7.69 

The following: table shows the steady state values of the transfer 
function matrix before and after optimization, (normalized by 
columns with diagonal elements being 100) 


Before Optimization After Optimization 



100. 

2.7 

6.7 

100. 

1.14 

9.9 

^ «5° 

19. 

100. 

8.8 

14.8 

100. 

4.4 


1.64 

2.1 

100. 

3.5 

0.69 

100. 


100. 

0.007 

0 

100. 

2.8 

3.1 

:^-10° 

0.001 

100. 

0 

6.1 

100. 

3.9 


0.007 

0.001 

lOQ 

5.6 

2.3 

100. 


100. 

9.7 

7.5 

100. 

2.0 

4.3 

^-15® 

3.7 

100. 

3.7 

9.4 

100. 

0.5 


2.6 

2.8 

100. 

8.02 

1 

5.0 

100. 


Time histories for the decoupled system with optimized (F,K) 
are given in the next four pages. It can be seen that there is an 
improvement in the steady-state coupling for the (2,1) element, 
which has been reduced from 19 to 14.8 for the'< “5° case. For other 
elements such as the (3,1) element in the o(* 15^ case, the steady 
state value has been increased about threo times. 

Based on these results, it seems that the decoupling matrices 
(F, K) has to be finely adjusted for the parameter variations In 
(A,B, C). Although we are able to find (F,K) through optimization 
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procedure so that some couplings are reduced, however, It seems to 
be very difficult to achieve significant reduc: lions for all couplings. 
Therefore, un alternative approach will be adc/ptod. Details are 
given in Section IV. 
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IV. Future research directions 

The following derivation is based mainly on C3p (see Chapter 

10 ). 

Let the system plant transfer function be denoted by P(s), 
which has m inputs and m outputs. Let E(s) and G(s) both be square 
rational matrices which are the compensators to be designed. 



Due to system parameter variations, the plant transfer function 
may become P'(s). Let L PG and L* ^ P'G. The over-all transfer 
function T - (I + L)“^ LE. It is easy to see that 

- E'^l'^ (I + L). 

E - L"^ (I + L) T. 

Let denote the over-all transfer function when the plant 
changes to P* , tT.e., 

T' - ( I + L')“^ L'E 

- E’^L'’^ (I + L') 




-1 



E"V"^ (I + L') 


-e”^lT^ (I + l) 


Hence 
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- e"^ (L'“^ - iT^). 

TT'“^ - I - TE"' (L'“^ - L"^) 

T - T* - TE“^ (L'“^ - L“^)T' 

Let iT ^ T' - T, then from the above equation and (5), 

AT - T {l"^(I + L)t} ^ (L“^ - L'“^)T' 

- (I + L)"^ (I - LL'“^)T'. 

- (I + L)”^ (I - PP'“^)T'. 

Suppose AT is relatively small, t.e,, T « T', then the above 
equation can be further reduced to the following, 


AT « (I + L)“^ (I - PP'"^)T (6) 

This equation is very useful in the case that the nominal 
plant transfer function P is diagonal. It is natural to assume L 
is also diagonal, then G ^ P'^L is also diagonal. Since both P and 
P^ are known, the matrix W ^ (I - PP^"^)T is also known. Consider 
the simple case where m - 2, then equation (6) can be rewritten as 


Mu 

Mj2 

«21 

“22 


rr 1, 


0 


0 

1 

1 + ij 


i 

j 

j 

j 


*^11 ^12 
W21 W22 


Suppose that the variations of T(s) are required to satisfy 
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|At|^j(.jw)| < kij(w), where kij(w) are sone given bounds on the 
elements of At^jCjw), then ^ 


Wij(.1w) 

1 + li(.}w) 


< kij(w). 


or 


Wjj(jw) 

jTij^w) ^ 


1 + Ij(jw) 


(7) 


The above equation shows that in order to keep the variations 
of T(s) within certain bounds, (i.e., kjj(wJ),, then the loop 
transmission L ^ diag (Ij, l£) has to be above certain value. 

Our next task is to find G ^ diag (gj, g2> such that L * PG 
satisfies the above requirement. Furthermore, all zeros of (1 1^) 

must be stable, i.e., the over-all system must be stable. This 
will be done on a cut and try basis on Bode plot. Finally, E(s) 
is choosen to realize T(s) as close as possible. 

In summary, the given system (A, B, C) will first l>e decoupled 
by using Gilberts* decoupling program, the decoupled system trans- 
fer function is called P. Then we design compensators E and G to 
reduce the effect due to parameter variations of (A, B) . 

New results on robust servomechanisms [4] should be very useful 
in the design of compensators E and G. 
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An Example in Decentralized Control Systems 

by 

Shih-Ho Wang 

Abstract - This note provides an example of a system which cannot be stabilized 
by decentralized state feedback, even though the over-all system and each local 
subsystem are completely controllable. 

In [ l] , Aoki and Li have made some interesting studies on the controllability 
and stabilizability of decentralized dynamic systems. They have suggested the 
following result: 

Consider a linear dynamic system described by 


3^-. = A_, X + A-_ X + u, . 
2 211 22 2 22 2 


^11 ^12 


Assume that the pairs (A, B.J, (A*,, B,_)and,,. . 

11 11 22 22 \ |A, , A-- 

6 1 66j p 


0 B 


22. 


are completely controllable, then one can find a set of decentralized (local) state 


feedback u. = C. x^, (i = l, 2) such that the resulting system matrix 




1 


21 


A.,, + B,, C, ! 
22 22 2J 


is stable, (see Proposition 7 of [ l] ). 
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Department of Electrical Engineering, University of Maryland. College Park, Md. 


20742. 


In the following example, we consider a system with three local control 


stations and the pairs (A, B), (A.., B_), (i = 1. 2, 3) are all controllable. Yet 


the system cannot be stabilized by local state feedback. 



, 0 0 

i 

i 1 0 

( 

I 

j 0 0 




0 

i 

10 

I 


; 0 

' 0 
I 









where is the control input at station i, (i = 1, 2, 3). The set of decentralized 
local state feedback is of the form 


u, 

1 


•'izl 



(i= 1,2,3). 


Applying the above set of feedback, the over-all system matrix is as 
follows: 



It can be easily seen that the first, second, third and fifth rows of the 
above matrix are linearly dependent for arbitrary other words, the 

above system has a fixed mode [2] at the origin. Hence the closed-loop system 
cannot be asymptotically stable. 
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An Approach to the Design of Decoupled 
Systems in the Presence of Large Pareuneter Variations* 

Shih-Ho Wang 
University of Maryland 
Department of Electrical Engineering 
College Park* Maryland 20742 

Peter Dorato (Tels 505/277-2436) 

Department of Electrical Engineering 
and Computer Science 
University of New Mexico 
Albuquerque, New Mexico 87131 

Summary 

In this study we outline a design approach -for developing 
control systems in the presence of large parameter variations. 
The decoupling problem is of practical interest in flight 
control problems, especially in the landing mode where it is 
desireable to independently control certain output variables, e.g. 
flight path angle, foward speed, etc. A common source of 
parameter variations in flight control problems is the result 
of linearization about a number of different operating points. 

It is assumed in this study that parameter variations may be 
represented by a discrete finite set of possible values, cor- 
responding to different flight conditions and that the complete 
system state is not available. 

The theory of decoupling with no parameter variations is 
well developed. See, for example, references [1] and [2]. 

Some work has been done by Fabian and Wonham, [3], on decoupling 


♦This research has been supported by NASA Langley Research 
Center under grant NSG 1213 at The University of Colorado. 


wits small parameter variations. However, very little vrork 
has been done on decouplin9 with large parameter variations. 

In reference (4) , Dora to, Wang, and Asher develop conditions 
for perfect decoupling with a discrete finite parameter set. 
However, in most practical cases the conditions for perfect 
decoupling are not satisfied. Thus in general, with parameter 
variations it is necessary to minimize some measure of inter- 
action in any decoupling design. 

A classical approach to the problem posed here is to 
use large loop gain aroung the plant to reduce the matrix 
transfer function to a good approximation of the unit matrix 
and then precede the feedback system with a known diagonal pre- 
compensator. An outline of this approach may be found in 
chapter ten of reference 15]. This approach has considerable 
merit, however, it is a trial and error approach which in 
addition requires a large number of graphical displays. 

The design approach suggested here involves the following 
steps. 

1. A full order observer, [6], is designed, based on a set of 
nominal parameter values, to generate a state estimated, x. 

2. A feedback control law of the form u = Fv + Kx is selected 
with design parameters consisting of all or selected 
entries of F and K. 

3. An interaction function is minimized with respect to F 
and K, subject to constraints on the diagonal elements of 
the closed-loop transfer function. The minimization problem 
is then solved via a nonlinear programming algorithim. 

The interaction function is selected to be of the form 

2 

F z r. |h..(ju))j 

a i,j ui 

i/i 



Where o denotes the parameter set* u) the frequency variable » 
and the i-j entry of the closed- loop transfer function. 

The above procedure involves a considerable amount of 
computation for any practical problem. However, all the steps 
are based on well documented algorithim which have been 
computerized. Some preliminary computational results indicate 
that the interaction measure can be reduced compared to a 
design based entirely on nominal values. 


References 


1. A. S. Morse and W. M. Wonham, "Status of Noninteracting 
Control", IEEE Trans* Automatic Control , vol. AC-16, 

no . 6 , December , 1971, pp. 568-581. 

2. E. G. Gilbert and J. R. Pivnichny, "A Computer Program for 
the Synthesis of Decoupled Multivariable Feedback Systems , 
IEEE Trans. Automatic Control , vol. AC-14, Dec. 1969, 

pp. 652-659. 

3. E. Fabian and W. M. Wonham, "Decoupling and Data Sensitivity", 
IEEE Trans. Autom atic Control , vol. AC-20, June 1975, 

pp. 338-344. 

4. P. Dorato, S-H Wang, and R. Asher, "Perfect Decoupling ^ 
of Linear Systems with Discrete Parameter Uncertainties , 

IEEE Trans. Automatic Control , vol. AC-22, no. 3, June 1977, 
pp. 498-500. 

5. I. M. Horowitz, "Synthesis of Feedback Systems", Academic 
Press, New York, 1963. 

6. D. G. Luenberger, "An Introduction to Observers", IEEE 

Trans. Automatic Control , vol. AC-16, no. 6, December 1971, 
pp. 596-602 . ~ 


491 


Itn TIUNlACnONI OH automatic CONTIOU lUHt I9T7 


QVAXJJ^ 


Perfect Decoupling of Linear System« with 
Discrete Parameter Uncertainties 

PETER DORATO, SHIH-HO WANG, and ROBERT ASHER 

ittefrarf — A flnplc dniga procedure b preicotcd, based ou GUben’s 
dccoupHag paramelcn (tH^l for perfect scaler d<coupliii| ria stale feed- 
back of a Uaetr lltue-lBvartaoi syslen with paraaieler uacenalailcs de- 
flacd over a discrete set of possiMe rahiet. 


Mtaiitchp* rtctivtd tuM 21, ISTS. Thu work wu lupporud ia pan by iht Naiwaal 
Atroaauuct and Spaca Adauniauaiioa tuidar Oraal NSO 1213. 

f. Doraw la with Iha Dtpartmanl of ElactrKal Eaftaatnai aod Compular Sctaaca, 
Uaivanily of Nto Maiico, AlbuAuarqua. NM ITIII. 

S.*H. Waag u anih Uia Dapanataai of Elacincal Eagiatanag. Uaivanity of Colorado. 
Colorado Spnaga CO lOtOt 

a. Aahar ia wilt Iha Dapartaiaal of Elacincal Eagiaaanag, Tatu Tacb Uaivarwiy, 
Lubbock. TX mat. 

Copyrighi ® 1977 by The Institute of Electficsl and Electronics Engineers, Inc. 
Printed in U.S.A. Annais No. 706AC048 



TECHNICAL NOTES AND COMESPONDENCE 


ORIGINAL PAG?: lb 

Or' Pi^'>K QUALr''V 


499 


I. Introduction 

The design problem it lo determine a fiaed slate feedback control law 
of the form 

u(/)-fjt(/)+Gt'(/) (I) 

which decouples the linear dynamical system, with input i (r) and output 

XO. 

I >(i)-C^.T(r) 

for all admissible system parameters The index; ranges over 

a finite set of integers J. and represents a discrete uncertainly set for the 
system matrices A. B. and C. This type of uncertainly appears naturally 
in systems which are linearized about a finite number of operating 
poinis. such as flight control systems and process control systems. In 
addition, this discrete model may be used as an approximation lo 
problems with continuous parameter uncertainties. The literature on 
decoupling of systems with no uncertainty is rather extensive. See, for 
example ll). (21. (b). |8]. and (9|. The literature on decoupling with 
parameter uncertainties is much more limited. In (31 the problem of 
decoupling in the presence of continuous parameter uncertainties is 
considered from the geometric point of view. Necessary and sufficient 
conditions are presented for the existence of a fixed state feedback law 
which preserses input-output properties of a system, including decou- 
pling and stability, for sufficiently small parameter variations. The 
present study differs from (3] in the following respects. 

1) The uncertainty set is discrete rather than continuous. 

2) No restrictions are placed on the size of the parameter variations. 

3) The diagonal elements of the decoupled system transfer matrix are 
not required to remain invariant or stable with respect to parameter 
variations. 

4) The output matrix C also may be variable. 

In (3| the invariance with respect to parameter variations of the 
mpui-outpui properties is referred to as data insensitivity The relaxation 
of the conditions of data sensitiviiy and stability for the diagonal 
elements in the closed-loop system permit the inclusion of a wider class 
of problems. Problems of invariance and stability of the diagonal ele- 
ments can be treated, if necessary, by single-loop techniques such as in 
(4). However, even with relaxed conditions on diagonal elements many 
practical systems cannot be perfectly decoupled, as illustrated in Exam- 
ple 2 in Section III 

In Section II conditions are given for perfect decoupling. These 
conditions lead di.ecily to the ■'equired control matrices F and G. 


II. L)E( Ul PLt^C Dt SKiS 

As shown by Gilbert (l|. any control maitices F and g w hich decouple 
the system (2). fur fixed >. must be of the form 

F--D 'A»*t 2 (0,,-ir,.)^ • V (3) 

m 

G- 2 \C, 

<•1 

where D. A*. p,. ir,*. 7,*, A, 4. and G, are parameters fixed by 

iAj.B^.C^) We refer to these parameters as the oilbert decoupling 
parameters. The parameters A,, a,^, and p,^ are free parameters which 
may lake on arbitrary real values. To stress the dependence of the 
control matrices and Gilbe'i decoupling parameters on (A^.B^.C^} we 
denote F as Flj), J,t as etc As further shown in (l|. the free 

parameters o,^ determine the poles of the closed-loop diagonal elements 
and A, determine the gain factors for the diagonal elements. A detailed 
program for the numerical computation of the Gilbert decoupling 
parameters is described in (2). We can now state the condition and 
design procedure for decoupling, given any index set 
The system (2) can be decoupled for all j belonging to J if and only if 
there exists a set of free parameters A,(/), and p^(j) for every y, 
such that 


FijJ (5) 

C(y,)-G(yj) G(y^). (6) 

If a suitable set of free parameters exists, the required control matnees 
are given by the common value of F and of G in (S) and (6). If the same 
free parameters can be used in each f(y) and G(y) to satisfy (S) and (6), 
then the system is data insensitive in addition to decoupled. Finally if the 
free parameters can be selected so that the poles of the diagonal 
elements all have negative real parts, the decoupled system also will be 
input-output stable. 

It should be noted that (S) and (6) represent a system of linear 
equations in the variables A,(y). o,i(y). and p.^fy). 


III. Examples 

Example I. Consider the system 

x^A^x + Bu, y^Cx, 7- (1.2.3) 


^.-(1 ;!)■ :;)■ ^>-(1 :)) 

*-(i !)■ ^-(i -1) 

The control matrices for this example are 

/"’■(-I ?) 

?) 

») 

The following values of free parameters satisfy (5) and (6). 

o,,(l)--l. 0||(2)--l. o,,(3)-0 
Oj,(l)--l. 02 ,( 2 )>O. Oj|(3)"-I 

A, ( I )- A, (21- A, (31-1. Aj(l)-Aj(2)-Aj(3)- I 
The resulting control matrices are 

;)■ =-(i ?) 

In this problem the diagonal elements are given by 

\(;) 

j-o,i(y) 


Note that the closed-loop poles of the diagonal elements vary with the 
parameter variations The condition for data insensitivity required in (3| 
IS not met by this example However, n can be perfectly decoupled for 
y-1.2.3. 

Example 2 The numerical values in this example are obtained from 
the lineaiized longitudinal dynamics of a STOL aircraft (7|. The uncer- 
tainty correspond! to two possible operating conditions, one correspond- 
ing to an angle of atuck equal to 3* and the other to an angle of atuck 
equal to 10*. For convenience we select (3, 10). 
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y-Cx 



0 

1 

0 

0 

0 

-1.285 

-0 6725 

0.209 

0 

1 

-0.403 

-0.570 

-0.286 

0 

0.1105 

-0.0785 

0 

0 

0 


0 

-2.315 

0.127 


0 

-0.0646 

-0.190 


0.1047 

-0.00383 

-0.1247 


0 

1 

0 

0 

0 

-1.23 

-0.520 

0.225 

0 

1 

-0.368 

-0.640 

-0.3195 

0 

0.157 

-0.1018 


I 0.12196 0 0 0\ 

•«-*u(IO)| -0.41003 0 0 Ol 
\ 0.I6I9I 0 0 0/ 

/O 0 0 9.5SII\ 

-»-«)i(IO) 0 0 0 0 I- 

\0 0 0 0 / 

Noic ihil the fret parimeiert ere miuing from the and entnet of 
the F{5) and f(IO) matnccs; hence no choice of free parameien can 
uiisfy f(S)> f (10). Thus, ihit lysiem cannot be perfectly decoupled for 
all admiuible parameter variations. The only reasoiuble decoupling 
design criterion in this case is to minimize some measure of interaction 
M in (S). 


0 

0 

0 

0 

-2.38 

0.1487 

0 

-0.0676 

-0.1712 

0.1047 

-0.01406 

-0.1190 

1 0 

-1 0\ 


1 0 

0 0 


0 0 

0 1/ 



The control matrices F(J) are in this case 


/•(5)-(o 

\o 


0.I8S03 

-0.30231 

0.17079 


-3.4599 
-0.36818 
- 1.99587 


- 2.79854 \ 

-0.06724 

-2.97714/ 


IV. Conclusions 

Conditions are presented for perfect decoupling of linear systems with 
ducretc parameter uncertainties. The conditions involve a test for the 
existence of a solution to a system of linear equations. An actual solution 
of the linear equations yields the fixed decoupling control law. While 
there exisu systems which can be perfectly decoupled in the presence of 
parameter variations, e.g.. Example I, many practical systems cannot be 
perfectly decoupled, e.g., Example 2. However, it is of some interest to 
know when perfect decoupling is possible. 
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/ 6.17229 0 - 6.17229 0\ 
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( 0.144 0 0 0\ 

-0.i9C9 0 0 0 

0.13291 0 0 0/ 

/O 0 0 9.5311 \ 

•t-e3i(5) 0 0 0 0 

\0 0 0 0 / 


and 


/ 3.05158 0.15862 -3.83436 

fdO)- 0 - 0.50436 - 0.34429 

Vo 0.19915 -2.01358 


- 3.23393 V 

-0.13568 

-3.68474/ 


/ 6.52689 0 - 6.52689 

■♦•o„( 10) 0.35616 0 - 0.35616 

\ 5.70049 0 -5.70049 

(0 0.12896 0 0\ 

•«-e„(IO) 0 - 0.41005 0 0 
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An Algorithm for Obtaining the Minimal Realization of a Linear 
Time- Invariant System and DctcrmininR if a System is 
Stabil liable- Detec table 

by 

* * t 

E.J. Davison , W. Gesing , S.ll. Wang 


ABSTRACT 

A definition of centralized fixed modes of a system, which is an adoption 
of the definition of fixed modes for a decentralized system [1], is made. It 
is shown that the centralized fixed modes of a system can be easily and effi- 
ciently calculated in terms of the eigenvalues of the system, and that the 
calculation of such fixed modes enables one to determine, in a numerically 
efficient may, whether a system is rpntrol lable , observable , stabilizable , 
detectable. An efficient algorithm for reducing a system to minimal realization 
form is then given. The algorithms have been effectively used on systems of 
order up to 100. 

1. INTRODUCTION 

Although the concepts of controllability, stabilizability , observability, 
detectability, etc. are well known for linear time-invariant systems, the 
actual determination whether a given system has one of the above properties 
is non trivial especially when the order n of the system is high (e.g. n>>10); 
this is because (i) most methods usually require the calculation of the rank 
of a certain matrix (which gets larger as the order of the system increases), 
and this always involves a decision, more or less arbitrary, as to what value 
to choose for zero (i.e. zero = 10 ^ or =10 • (it) often the rank of a 

large number (n) of matrices of order n must be found or the rank of a single 
matrix, (of order n) , but ill-conditioned, must be found. 

nxn n^m 

For example, it is well known that the pair (A, B) c R x r is 
controllable if: 
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( 1 ) 


rank (B, AB,.*.., a"‘‘b) = .. 

In this case, however, the determination whether (1) is true will be extremely 
difficult for large n (e.g. n»10) due to the ill-conditioning of the con- 
trollability matrix (B, A!., ..., A^'^B) (i.e. the columns of the controllability 
matrix tend to become linearly dependent for largo n ). Alternately, it is 
well known that the pair (A, B) is controllable if: 

rank (A-X^I, B) = n , i = l,2,...,n (2) 

where X., i = l,2,...,n are the eigenvalues o' A [2J. In this case, however, 
the rank^of a nx(n+ra) matrix must be determined n times, and this is non 
trivial to do for large n . It will be shown that the above problems can be 
eliminated by calculating a set of p scalars for the system (0 < p n) 

(called the centralized fixed modes of the system), which only requires the 
determination of the eigenvalues of two nxn matrices for the system. 

.Similarly, although the concept of the minimal realization of a linear, 
time-invariant system is well-known, the actual determination of a minimal 
realization for a linear system is non-trivial. especially when the order 
n of the system is high (e.g. n » 10) . This is because most existing algor- 
ithms, e.g. [3]- [8], involve a branching process, which depends on whether 
a certain number is zero, and subjective decisions must be made as to what 
value to choose for zero. Different realizations may thence obtained^ 
depending on what value to choose for zero (i.e. zero t 10 or t 10 )• 

It will be shown that the above problem can be eliminated if the centralized 

fixed modes of the system are known. 


2. DhVELOPMENT 

The followins definition is required in the development to follow; it 
is adopted from Wan.;; and Davison [I), who defined fixed modes of a decent rali-.ed 

system. 

Definition 

consider the triple (C.A.B) e « R"”" = then the set of 

centralized fixed modes of (C,A,B) , denoted by A(C.A.B) is defined as 

fol lows: 

A(C,A,B) « n o(A*BKC) 

KcR 

where o(A>BKC) denotes the list of eigenvalues of (,UBKC). 
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Remark 1 

The set of centralised modes of a given triple (C,A,Bj can be calculated 
as follows. Let k(i, j) denote the (i, element of K. Since K 

contains the null matrix, then A(C,A,B) c o(A) = . For each 

€ o(A), det (X^I - A - BKC) is a polynomial in k(i, j), i = l,2,...,m ; 

j » 1,2 r , and if dot (X^^I - A - BKC) is identically zero, then clearly 

X^ € A(C,A,B) . On the other hand, if .det (X^I - A - BKC) is not identically 

^ niXT* 

zero, then there exists K e R so that det (X^I - A - BKC) i 0 , and 
hence X^ ^ A(C,A,B) . This procedure determines the fixed modes of A(C,A,B). 

The following lemma is obtained by using identical arguments as used in 
the proof of theorems 1-3 of (9]. 

Lemma 1 

Given (C,A,B) t R^’*" x r"’'” x r"*'"’, then the class of matrices K e 
which does not result in the centralized fixed modes of (C,A,B) being equal to 
those eigenvalues of A+BKC which are common with the eigenvalues of A, is either 
empty or lies on a hypersurfacc [9] in the parameter space of K. 

The above lemma forms the basis of the following algorithm to deter- 
mine the centralized fixed modes of (C,A,B), 

Algorithm I: to Determine the Centralized Fixed Modes of (C,A,B) 

1. Find the eigenvalues of A . 

2. Select an arbitrary matrix K € R*** (by using a pseudo-random number 

generator, say) so that || Ajj = || BKCH . 

3. Find the eigenvalues of A+BKC . 

4. Then, the centralized fixed modes of (C,A,B) are contained in those 
eigenvalues of A+BKC which are common with the eigenvalues of A (if 
there are any). Moreover for "almost any" K chosen, the centralized 
fixed modes will be identically equal to tliose eigenvalues of A+BKC 
which are common with the eigenvalues of A , i.e. the class of matrices 
K which does not result in the fixed medes being identically equal to 
the eigenvalues of A which arc common with the eigenvalues of A+BKC 
is cither empty or lies on a hypersurfacc [9] in the parameter space of 
K . 

5. Repeat steps 2 to 4 using a different K if there is any doubt which 
eigenvalues of A arc the fixed modes of (C,A,B). 
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Proof 

The proof of statement 4 in the above algorithm follows directly from 
lemma 1 . 

Remark 2 

The only possible source of difficulty which may occur in the above 
algorithm is in making a decision as to what eigenvalues of A arc common 
with the eigenvalues of A+BKC . For example, it may occur that an eigen- 
value of A only differs from an eigenvalue of A+BKC in the 7th significant 
figure. Such cases can normally be decided upon by utilising the physical 
properties of the plant associated with the mathematical model (C,A,B). 

For example, if the plant data is only accurate to 3 significant figures, 
then in the above case, it can safely be assumed that the eigenvalue in 
question is a fixed mode. 

3. APPLICATION OF CENTRALIZED FIXED MODES 

An application of the above algorithm can immediately be made to deter- 
mine certain properties of (C,A,B) and to find the minimal realisation of 
(C,A,B) . The following results follow directly from the definition of 
centralized fixed modes: 

To Determine if (A, B) is Controllable 

(A, B) is controllable if and only if (I^, A, B) has no centralised 
fixed modes. 

To Determine if (A, B) is Stabilizable 

(A, B) is stabilizable if and only if (I^, A, B) has no centralized 
fixed modes which lie in the closed right half part of the complex plane. 

To Determine if (C, A) is Observable 

(C, A) is observable if and only if (C, A, 1^) has no centralized 
fixed modes. 

To Determine if (C, A) is Detectable 

(C, A) is detectable if and only if (C, A, I^) has no centralized 
fixed modes which lie in the closed right half part of the complex plane. 

To Determine if (C, A, B) is Controllable-Observable 

(C, A, B) is controllable and observable if and only if (C, A, B) has 
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no centralised fixed modes. 

To Determine if (C. A, R) is Stahi 1 isahle-Dctcctable 

(C, A, B) is stabilizable and detectable if and only if ((', A, B) has 
no centralized fixed modes which lie in the closed right half part of the 
complex plane. 

Remark 3 

The centralized fixed modes of correspond to 

the input-decoupling zeros, output-decoupling zeros of (C, A, B) respectively, 
defined by Rosenbrock [2]. 

To Determine the Minimal Realization of (C, A, B ) 

The following algorithm is obtained: 

Algorithm II 

1. Determine the centralized fixed modes of (C, A, B) ; let these modes be 
denoted by {Xj ,X 2 » • • • .^p> . 0 ^ p s n . For non-trivality, assume 

p < n . 

2. Apply a arbitrary output gain matrix K to the system (by using a 
pseudo-random number generator, say so that || a| 1 = || BKCH ), so that 
the controllable-observable modes of (C, A, B) , A = A+BKC are all 
distinct and disjoint from the fixed modes of (C, A, B) . This will 
be true for "almost all" gains K from [9]. 

3. Determine the eigenvalues of A . Let these eigenvalues be denoted by 

{Xj.X^,. .. ,Xp, * 

4. Determine the eigenvectors of A for the eigenvalues X^, i = p+1 ,p*2, . . . ,n . 
Let these eigenvectors be denoted by x^, i = *p-*-l ,p+2, . . . ,n respectively. 

5. Determine the eigenvectors of A’ for the eigenvalues X^, i = p+ 1 ,p+2 , . . . ,n 
Let these eigenvectors be denoted by y^^, i = p+1 ,p*2, . . . ,n respectively. 

6. Normalize the eigenvectors y^, i = p*l ,p*2, . . . ,n so that ~ 

i = p+1 ,p*2, . . . ,n . 

7. Then the minimal realization of (C, A, B) has order n-p and is given 
by (C, A, B) where: 
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^ 

A « 

f 5p.I ] 

yp*2 


(4) 


y’ 


n 


and the minimal realisation of (C, A, B) is given by (C, A-BKC, B) . 

The above algorithm produces a minimal realization which in general has 
complex elements appearing in (C, A, B) ; the following algorithm gives a 
minimal realization which has only real elements appearing in (C, A, B) . 


Algorithm III (to Determine Real Minimal Realization of (C, A, B) ) 

Assume that steps 1 to 5 in the above algorithm have been carried out, 
and assume that the ordering of the eigenvalues of A has been carried out so 
that: 

« 

^^p+l’^p*2’“’’^n^ " ^^p+r^p+2’ ” ' ’^p+t' ^p+t+r^p+t+r^p+t+2’^p+t+2” * 

...,X ^ ,X ^ } (5) 

* p+t+q p+t+q 


are 


where Xp^^ Xp^^ are all real eigenvalues, Xp^^^j ,Xp^^^^, ... ,Xp^^^^ 

all complex eigenvalues and t+2q = n-p . Steps 6, 7 of the previous algorithm 
now become: 

* A ^i 

6’ For i = p+1, p+2, ..., p+t, let y^ = . 

1 Xi^i 

For i = p+t*l, p+t+2, ..., p*t*q, let y* = Rc(y^) - Inj(y^)] + 

■* a ♦b a +b“ 

j[- — ^Re(y.) - I">(y )]* '^*'cre a = Re'(yJRe(x ) ♦ Im' (y . ) Im(x. ) , 

a +b ^ a *b ^ 


b = Re' (y^)Im(x^) - Im' (x^)Re(x^) . 
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I.ct r, . f, /n-p)Mn-p) follows: 

r, = {x ,,x -,...,x /2Rc(x . ,),/2lm(x . ,),/2Ilcfx . _),/2Ini(.\ ) 

1 ‘ p-^1 p^2 p*t p+t^l p+t^r p^t+2 P^t*2 


(7) 


A = block diag 


V2’ 


..X 


p+t’ 


Rc(Xp^^^l). Im(Xp^^^j) 




Re(X , ), Im(X ^ ) 

p+t+q P+t+q 


Iin() , ). Re(X ^ ) 

p+t+q P+t+q 


( 8 ) 


Then the minimal realization of (C, A, B) has order n-p and is given by (C, A, B) 
where C = Cr^, B = T2B, and the minimal realization of (C, A, B) is given by 
(C, A-BKC, 6). 

Proof 

The proof of statements 7, 7' in the above algorithms follow directly by 
applying a Jordan decomposition to (C, A, B) and using the definition of fixed 
modes of (C, A, B) . The details arc omitted. 

Discussion 

In step 7' of algorithm III, the matrices Fj, F^. A arc all real, so that 
the minimal realization (C, A, B) obtained contains only real elements. 

In steps 3-5 of algorithm 1 1 and 1 1 1, the calculation of the eigenvalues 
and eigenvectors of a matrix with only distinct eigenvalues is required. It 
should be noted that many extremely effective algorithms exist to do this 
(c.g. the QR family of algorithms) and are readily available (e.g. the RISPACK 
set of subroutines). It should also bo noted that the computation time required 
to carry out all steps of algorithm 1 1 1 only varies as n^ (the computation time 
required to carry out steps 3-5 of algorithm 1 1 , 1 1 1 only varying as (n-p) .) 
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lixam])]c 1 

It is desired to find a state space minimal realization of a boiler furnace 
model 1 10) whose transfer function matrix is Riven by: 


G(s) = 


1 

0.7 

0.3 

0.2 

1^4s 

l^Ss 

U5s 

U5s 

0.6 

1 

0.4 

0.3 

l-»5s 

l+4s 

l+5s 

l+.Ss 

0.35 

0.4 

1.0 

0.6 

l^Ss 

l+5s 

l+4s 

l+5s 

0.2 

0.3 

0.7 

1.0 

l+5s 

l+5s 

l+5s 

l+4s 


(9) 


In this case, a 16 order nonminimal rtate space realization of (9) is directly 
obtained as: 


Hj = -0.25 Hj + 0.25Uj 
= -0.20H2 ♦ 0.14u2 

Hj = -0.20nj 0.06Uj 

n* = -0.20n. + 0.04u 
4 4 4 

Tig = -0.20nj ♦ 0.12Uj 
n, = -0.25n, ♦ 0.25u- 

O Q Z 

= -0.20Hy ♦ O.OSUg 
fig = -0.20ng ♦ 0.06u^ 


Tig = -0.20ng ♦ 0.07Uj 

^10 ' 

n,, = -0.25n,- + 0.25u, 
11 11 a 


'12 


-0.20nj2 ♦ 


^13 * 

^14 “ ♦ O.O 6 U 2 

fijg = -0.20 t1j5 ♦ 

"j6 ■ -“•^‘”’16 * 


( 10 ) 


>'l ■ "1 * "2 * "3 * "4 

>'2 ■ "5 * "6 * "7 * ''8 

)'3 ■ "9 * "10 * ”11 * ''12 

>'4 = "l3 * "l4 * "is * "16 

and on applying the feedback u = Ky to (10), where K is given by (arbitrarily 
chosen) : 


-.368462 

.255606 

-.041349 

.032768 

-.28104 

-.452955 

.178865 

.179297 

.434693 

-.116497 

.019417 

. 330966 

-.465427 

-.446538 

.029701 

.17115 
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the following open loop anJ closed loop eigenvalues of (10) were obtained 


(see Table 1): 


Table 1 Open and Closed Loop Eigenvalues of (10) with u = Ky 


nicenvalucs of A 

nigcnvalucs of A+BKC 

-0.200000 

*-0.200000 

-0.200000 

*-0.200000 

-0.200000 

*-0.200000 

-0.200000 

*-0.200000 

-0.200000 

*-0.200000 

-0.200000 

*-0.200000 

-0.200000 

*-0.2000o0 

-0.200000 

*-0.200000 

-0.200000 

-0.343009 ± jO. 0571058 

-0.200000 

-0.236782 ± jO. 0630918 

-0.200000 

-0.225945 

-0.200000 

-0.164543 

-Q. 250000 

-0.170378 

-0.250000 

-0.201203 . 

-0.250000 


-0.250000 



*Fixed modes of (10) 


which implies that the fixed modes of (C.A.B) for (10) are given by (-0. 2, -0.2. -0.2, 
-0.2, -0.2, -0.2, -0.2, -0.2), and thus that the minimal realization of (10) 

has order 8. 

Remark 4 

It should bo noted that the eigenvalues of the open loop and closed loop 
system corresponding to the fixed modes agree with each other to 15 signi- 
ficant figures (using double precision arithmetic on an IBM 370 digital 
computer) . 

On applying algorithm III to the system (10), the minimal realization 
of (C, A, B) is directly obtained and is given (truncated to 6 significant 
figures) in Table 2. The minimal realization of (C, A, B) is then imniediate- 
ly obtained and is given in Table 3 (truncated to 6 significant figures). 


Remark 5 

It should bo noted that the transfer function representation obtained 
for (C, A, B) agrees with (9) to 15 significant figures. 
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Uxumplc 2 

It is desired to determine if the system (10) is control In this 

* * 

.case, on applying the feedback u = K x to (10), where K is given by (arlii- 
trarily chosen): 


-0.368 

-0.492 

-0.452 

-0.063 

0.256 

-0.116 

0.237 

0.267 

-0.041 

-0.433 

-0.171 

-0.022 

0.033 

-0.062 

0.133 

-0.262 

-0.261 

0.167 

0.257 

0.225 

-0.452 

0.089 

0.492 

-0.140 

0.179 

0.431 

-0.134 

-0.333 

0.180 

0.347 

-0.252 

-0.013 

0.435 

0.027 

0.483 

0.398 

-0.116 

-0.408 

0.223 

0.410 

0.020 

0.154 

0.254 

0.439 

0.331 

0.064 

0.152 

0.40S 

-0.465 

0.202 

-0.427 

0.005 

-0.446 

0.411 

0.132 

0.017 

0.030 

0.263 

0.385 

0.160 

0.172 

-0.237 

-0.227 

0.487 


'he following open loop and closed loop eigenvalues of (10) were obtained 
(Table 4): 

* 

Table 4 Open and Closed Loop Eigenvalues of (10) with u = K x 


Eigenvalues of A 

* 

Eigenvalues of A-»BK 

-0.200000 

•-0. 200000 

-0.200000 

*-0.200000 

-0.200000 

*-0.200000 

-0.200000 

*-0.200000 

-0.200000 

*-0.200000 

-0.200000 

*-0.200000 

-0.200000 

*-0.200000 

-0.200000 

*-0.200000 

-0.200000 

0.042654 

-0.200000 

-0.426849 

-0.200000 

-0.227648 ± jO. 072282 

-0.200000 

0.166892 

-0.250000 

-0.220722 ± jO. 004844 

-0.250000 

-0.212503 

-0.250000 


-0.250000 



*l'ixcd modes of (10) 


which implies that the fixed modes of (I^,A,B) for (lu) are given by {-0.2, -0.2, 
-0.2, -0.2, -0.2, -0.2, -0.2, -0.2), and thus that the system (10) is not 
controllable. 
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5. CONCIJISIONS 

A definition of centralized fixed modes of a linear, time- invariant 
multivariable system is made, and an algorithm (algorithm I) for calculat- 
ing the fixed modes of a system is given. The calculation of the fixed 
modes only requires the calculation of the eigenvalues of two n><n matrices 
for a n^*' order system. It is then shown that the determination of whether 
a system is controllable, observable, stabilizahlc, detectable etc. can be 
easily carried out in a numerically efficient way by calculating the cen- 
tralized fixed modes of the system. A direct a’^plication of fixed modes 
of a system is then made to obtain an extremely efficient algorithm 
(algorithms II, III) for obtaining the minimal realization of a linear 
time invariant system. The algorithms have been successfully used on 
systems up to 100^*' order. 
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2 Minimal Realitition of (C. A, B) Ohtiined For Exar.ple 


343009 

-.0571058 

0 

0571058 

-.343009 

0 

0 

0 

-.236782 

0 

0 

.0630918 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 0 

0 0 

-.06309IS *0 

-.236782 0 

0 -.225945 

0 0 

0 0 

0 0 


0 

0 

0 

0 

0 

-.164543 

0 

0 


0 

0 

0 

0 

0 

0 

-.170378 

0 


0 

0 

0 

0 

0 

0 

0 

201203 


’ .27667 

.550161 

-.0783017 

-.0686386 





.734103 

.102412 

-.0332067 

-.0851425 





-.0946539 

-.352436 

.263639 

.57111 





.212113 

.463087 

.124411 

-.223794 





-.0280378 

-.058489 

-.00390984 

.0748103 





-.36691 

.037034 

.433436 

.350038 





-.413039 

.124359 

.215827 

.178217 





.0268964 

.059385 

.00923714 

-.0262933 





.190429 

.266941 

.12383 

.123835 

.0553281 

.123578 

-.0813695 

.00434014 

.431295 

-.0556493 

.17273 

.20436 

-.012404 

.201877 

-.215097 

-.0209016 

-.180345 

.0613934 

.36793 

.665509 

.037844 

.181602 

-.0917767 

-.0614214 

.436308 

-.0539819 

.347657 

-.0890121 

-.199607 

.386107 

-.355563 

-.00109043 



Table 3 

Minimal Realization 

Of (C. A. 

B) Obtained For Examolc 


0.0038101 

O.OOJ35373S 

3.0111167 

0.0736151 

0.0177717 

*3.0111)717 

3.91)67787 


*3.757011 

0.03716711 

0.3301366 

0.071)753 

*3.03006787 

0.0716:08 

0.1361916) 


0. OOKOOMil 

*0.736376 

0.006 37603 

0.0761016 

*3.0086)303 

0.01517105 

*0,13531)37 

0.0!2;S37 

O.oi07l‘*7 

0.3333J7 

*3. 756611, 

0.0171633 

*3.9301656 

9. 1701'. 0 

0,11)16)11 

0.noc7<ikC7 

0.00?k1873 

0.13'’.6Cg33 

1. 03170-' J1 

*0.735)63 

1.09551716 

*9.996C1)73 

*3, 1'*1''766 18 

0. 4S 

*3.3ke6110 

*3.0670661 

3.317171 

0.1365'')! 

*0.737)73 

1. 317701) 

*0.03651677 

o.csii'jm 

*O.OltC 3503 

*0.1763536 

1. 0151756 

3. 00305766 

*0.11)8117 

*9.17)667 

*3.11 1110' 1 


O.OQ333»6« 

*0.33300311 

*0,10101736 

0.00711)3 

*0.907)677) 

0.001551)6 

*0.791705 

0,J7^^T 

9.5S0ir.i 

*0.0703017 

*0.96iC146 




• 

7.73<i1<»3 

0.107kl3 

*0.0113007 

*3.3351675 




. 

*3.0 

*0.357*136 

0.781653 

0.57111 





0.7i;i!3 

0.6C50C7 

3.136611 

*0.771736 





'3.0J7757* 

*<l.n5Mti 

*0.03510 306 

0.0760135 





~3.3bb71 

0.05733H 

0.611616 

0.357010 



. 


"3.«i>39»0 

0, 136351 

0.715037 

0.17i:i7 



* 

• 

0. .17bk«(ik 

0.051335 

0.00375716 

*0.3767331 





r 9.1331171 

0.768161 

0.17305 

0.171035 

0.055)761 

0.17)570 

*0.001)615 

1.006 36116*1 

|^3.ii3l3‘iS 

*1.0556613 

0.17773 

0,70616 

*'0.113606 

0.701077 

*3.715117 

*1.070)010 1 

|*3.1l)J3<it 

0. 0613136 

0.3..711 

0.6..i6 11 

0.117066 

9.101617 

*3.0117767 

'9,1':167!6 I 

L 6.*3f.339 

*3.0533011 

0.367857 

*0.04)3171 

*0.113607 

0.106107 

*0.15556) • 

9.9110336lJ 
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An Algorithm for the Calculation of Transinissiop Zeros of 
the System (C,A.B,P) Using High Gain (Output Feedback 

by 

• t 

E.J. Davison . S.H. Wang 

ABSTRACT 

A new algorithm, which is an extension of algorithm II given in [1], is 
given to determine the transmission zeros of the system x = Ax ♦ Bu, y = Cx ♦ 
denoted by (C,A,B,D) . The algorithm is based on the observation that for 
nondegenerate (C,A,B,D) systems, the transmission zeros of (C,A,B,D) are 
contained in the finite eigenvalues of the closed loop system matrix 

{a ♦ BK(-^ - DK)”^C} , where K is any arbitral/ matrix of full rank, > * R 
and p -► » . 


1. INTRODUCTION 


This paper is concerned about the efficient calculation of the trans- 
mission zeros of the following linear, time-invariant system denoted by 


(C,A,B,D): 


X = Ax ♦ Bu 
y • Cx ♦ [)u 


( 1 ) 


where x e r” is the state, u c r"* is the input and y f. R*^ is the output. 
It is not necessarily assumed that (1) is controllable or observable nor t!*.at 
rank B = m , rank C » r necessarily. 

In [l]-[3] the transmission zeros of (1) were defined to be the set 
of complex numbers ). which satisfy the following inequality: 


rank 


A- XI 
. C 


! < n ♦ nin(r,n) 

r>J 


(2j 


in particular, 


the transmission zeros arc the zeros (niultinlicitics i iiciuJ ctl. 
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of the greatest common divisor of all (n ♦ min(r ,mj) ><(n ♦ min(r,ni)) minors 
A-XI b] 

of . Two algorithms for calculating the transmission zeros of 

c dJ 

(1) were given in [1]. The purpose of this paper is to give a new algorithm 
for finding the transmission zeros of (1), which is an extension of algorithm 
II given in [1]. 

2. DEVELOPMENT 

The following definition is required in the development to follow: 
Definition [1] 

(C,A,B,D) is degenerate if for any specified n+1 distinr/ real scalars 
i = 1 ,2 , . . , ,n-*-l , the following is true: 

A-X*I , b] 

rank < n ♦ min(r,m) , i = l,2,...,n-^l (3) 

C , D 

A system which is not degenerate is called a non-degenerate system . A degenerate 
system has transmission zeros which include the whole complex plane. 

The following lemmas follow by observation: 

Lemma 1 

Assume that (A-XI) is nonsingular; then 

A-aI , , 

rank = n rank [D - C(A-XI) B] (4) 

. C , D 

Lezi ma 2 

Assume there exists a scalar X e I so that: 

A-XI b' 

rank = n + mini.r,ni) ; 

C dJ 

then (C,A,B,nj is non-degentTato. 

Lemma 5 

Given D c p/*"*, K c witn lanl. K -- r.infr.m) , then there exists 

I 

a constant a* > 0 so that (— - I'i:’! !'• ncaa i iniular Vp > p* . 

0 

* Ii,should be noted that an intere.->ung algorithm i.liidi is quite different from 
the one presented here has recently beuui suggested in [o]. 
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The following lenuna is required in the development to follow: 
Definition 

Given the polynomial equation in \: 
p-1 q _ 

C(XP ♦ I c A^) + e( [ c (c)X^) = 0 , q s p 

1 _ r\ * ' z ^r\ ^ 


( 5 ) 


i=0 


i=0 


with real coefficients c,c^,c^ with c / 0 such that: 

Cq(e) / 0, |c^(e)| < M, i = l,2,...,q, Vc €[0,»), let the roots of (5) be 
denoted by A., j l,2,...,q. Let the roots of the polynomial equation in 

P-1 i 

c(xP ♦ [ c A ) = 0 

i=0 

be denoted by Oy j = l,2,...,p. 

Lemma 4 [1] 

For any (6, "6), 6 > 0, T > 0, there exists an £*(6,6) > 0 so that 
Ve c[0,e*), p roots of (5) suitably ordered, have the property that 
[a. - Oj 1 <6, j = 1,2, . . . ,p and q - p roots of (5) suitably ordered, have 
the property that | A ^ | >7, j = p+1 , p+2, . . . ,q. 

The following lemma follows by observation from the definition of 
transmission zeros given in (2) . 


Lemma 5 

m><r 

Assume that (C,A,B,D) is nondegenerate and that K £ R has 
rank K = min(r,m) . Then: 


(i) 

If 

m > r 

for "almost all" 

K , the system 

(C,A,BK,DK) 

is 

nondegenerate 

(2) 

If 

m < r 

for "almost all" 

K , the system 

(KC,A,B,KD) 

is 

nondegenerate 

i.c.. 

the 

class 

of K matrices in 

which the above 

statements 

are 

not true 


is either empty or lies on a hypcrsurface [4] in the parameter space of K . 
Lem ma 6 

Assume that (C,A,B,D) is nondegenerate and that K c R has rank 
K = min(r ,m) . Then; 

(1) If m > r, for "almost all X", tlu transmission rcros of (C,A,B,D) are 
contained in the transmission zeros of (C,A,BK,DK). 

(2) If r < m, for "almost all K", t!n‘ tiviismission zeros of (C.A.B.D) are 
contained in the transmission zeros of (KC,A,B,KD). 

i.c., the class of X matrices in wliich the above statements are not true is 
cither empty or lies on a hypersurface in the parameter space of K. 
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Proof 

The proof of this result follows from lemma 5 on using the definition 
of transmission zeros given by (2). 

3. MAIN RESULTS 

The following initial result is obtained: 

Theorem 1 

Assume that (C.A.B.D) is nondegenerate; then for "almost all" 

X e I . rank (D-C(A-XI)'^B] = min(r,m), i.e. the class of X which results 
in rank [D-C(A-XI)"^B] < min(r.m) or in (A-XI) being singular is either 
empty or lies on a hypersurface [2] in the parameter space X e C. 

Proof 

The proof follows directly from lemmas 1, 2 on using iclentical arguments 
as used in the proof of theorems 1 , 2 of [4], 

Definition 

Given the system (1), let o(p) denote an eigenvalue of 

M(p) = {A ♦ BK(— - DK)‘^C} 

P 

where K e with rank K = min(r,m) and p e R . Assume now that there 

exists a scalar X* e I with the property that for any c > 0 there exists 
p(€) >0 so that io(p)-X*| < e Vp > p(e) ; then X* is called a finite 
eigenvalue of M(p) . 

Remark 1 

From lemma 3, the matrix M(p) in the above definition is well defined 
if p is sufficiently large. 

The following main result is now obtained: 

Theorem 2 (High Gain Output Feedback Theorem) 

Given the nondcgcncratc system (C,A,r>,D), let K c R be an arbitrary 
matrix with rank K = min(r,m). Then if m - r, the finite eigenvalues 

of {A ♦ BK(— - DK)”^Cl as p -> «> coincide with the transmission acros 
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of (C,A,B,D): if m r then for "almost all" K, the transmission zeros 

I 

of (C,A,B,D) are contained in the finite eigenvalues of (A ♦ - DK) C}, 

i.e. the class of K matrices in which the above statement is not true is 
either empty or lies on a hypcrsurfurce in the parameter space of K. 


Remark 2 

This theorem is an extension of theorem 8 of [1] and is a generalization 
of the classical result for single-input, single-output systems, in which it 
is well known that for high gain output feedback, the poles of the closed loop 
system approach the zeros of the system. 


I roof 


.The proof of this result follows very closely the proof of theorem 8 

in [1]. In particular following [1], assume with no loss of generality 

(theorem 7 [1]), that (A,B) is controllable and that (A,B) is in Brunosky's 

canonical form. Let n. , i = l,2,...,m be the controllability indices [5] 

5 ^ A 

of (A,B) with I n. = n. Let C = (C^ ,C 2 , . • • ,CJ be the output matrix 
i=l 


(not necessarily unique) which results from transforming fA,B) into this 

A ”i r*! 

canonical form and let D ■ (Dj ,D 2 » • • • where £ R and £ R 

Then the transmission zeros of (C,A,B,D) are given by the set {A} which 
satisfy the condition 


rank 


A->.I , B 
C . D 


< n ♦ min(r,m) 


which after substitution and simplification may be written as 


rank r(X) < min(r,m) 


( 6 ) 


where 


"(Cj.Dj) 

' 1 

, (L-^.n,) 

/ ^ 
1 

, . . . , (C_ * 

m m 

’ 1 



A 


A 


A 



• 

* 


• 


• 





• 


. 





A ", 


n 

A 

k J 



(7) 


j Assume now that m = r and that p i.s sufficiently large so that 
(“ - DK) is nonsingular (lemma 5;. Con.^iiler the eigenvalues of the matrix 

^ I 

M(p) = (A ♦ BK(— - DK) they arc given by the solution for \ of the 

P 

following equation: 

det (M(p) - AI) * 0 (8) 


. 6 


which can be simplified to the following equation (see Appendix 1); 

n. n_ n 

det {p K r(X) - diag(X )} = 0 (9) 

and on expanding the determinant of the left hand side of (9), the following 
is obtained: 

_ n-1 . 

p® det (K r(X)l ♦ {X ♦ I c.(p)X^) = 0 (10) 

i=0 

where the coefficients c^(p) have the property that 

lim fiiPl = 0 , i = 0,1 n-1. (11) 

p-*<» m 

P . * 

Since K is nonsingular, (10) may be written as; 

det (r(X)J ♦ {x" + I c (p)X^) = 0. (12) 

p®det K i=0 

Now since it has been assumed that (C,A,B,D) is nondegenerate, then 
det [r(X)] is a nonzero polynomial in X and hence from lemma 4, it follows 
that the finite eigenvalues of the matrix M(p) ^^e given by the solution 
of the equation; 

det [r(X)l = 0 (13) 

for X. This means, from (6), that the finite eigenvalues of the matrix 
M(p) are equal to the transmission zeros of (C,A,B,D) for the case m = r. 

The proof of tlio case m r follows directly from this result on using 

leiiimas 5,6. 

4. AN ALGORITILM FOR FINDING TR.ANSMISS ION ZEROS OF A SYSTEM 

Theorems 1, 2 form the basis of a new algorithm for finding transmission 
zeros of (1). The following algorithm is obtained from theorem 1. 

Algorithm to Determine if (G.A.B.P) is Non-degenerate 

f 

1. Choose an arbitrary real X* (by using a pseudo-random number generator 

say) so that 1| j || a|1 ) so that (A - .\*I_^) is nonsingular. 

This will be true for "almost all" X* chosen. 

2. Calculate rank (D - C(A-X*l)’^B) . If rank (U - C(A-XM) ^B) = min(r,m), 

the system (C,A,B,Dj is nondegencratc; if rank (D - C(A-X*I) B) < min(r,m), 
then (C.A.B.D) is "almost always" deecnerare . i.e. (C,A,B,P) can only be 
nondcguaiMc provided X* lies on a hypersur face in the parameter space 

X* c a . 
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3. If in doubt as to whether (C,A,B,D) is degenerate, repeat step 2 using 
a different real X* ; if rank (D - C(A-X*I) ^B) < min(r,m) for 
min(r,m)+ 1 distinct real values of X* , the system (C,A,B,D) is 
degenerate. 

The following algorithm follows from theorem 2. In this algorithm, it 
is assumed that (C,A,B,D) is nondegenerate, but it is not necessarily assumed 
that (C,A,B) is controllable or observable or that rank B = m, rank C = r 
or that m = r. 

Algorithm to Find Transmission Zeros of a Nondegenerate (C,A,B,D) System 

1. Choose an arbitrary matrix K c R™ ^ (by using a pseudo tcindom number 
generator say) so that rank K = min(r,m). 

A ^ -1 

2. Calculate the eigenvalues of the matrix M(p) = {A ♦ 

large p (e.g. p = 10^^ ) ; then if ra = r, the transmission zeros of 
(C,A,B,D) are equal to the finite eigenvalues of M(p), and if m f r, 
then for "almost all" K, the transmission zeros .of (C,A,B,D) are 
contained in the finite eigenvalues of M(p) • 

3. If in doubt as to what the finite eigenvalues are, repeat steps 1, 2 
with a different K and different p (e.g. p = 10^^). 

Remark 5 

In case the condition m = r is not satisfied, those finite eigen- 
values of M(p) which do not satisfy (2) should be discarded; the remaining 
finite eigenvalues arc then the transmission zeros of the system. 

Remark 4 

Note that it is not necessary to actually calculate the rank of matrices 
in the above algorithms, e.g. let .N = (D - C(A-X*1) 3) — then if r s m, 
rank N = r if and only if NN" is nonsingular, and if r > m, rank N = m 
if and only if N'N is nonsingular. 


* 

More precisely, 


p « 10 


15 



Here, it is assumed that double precis ior. 


arithmetic is used with this value of p . Tf 
used, a smaller value of p is suggested, say 
of rounding errors in the digital computer. 


.single precision arithmetic is 

g 

p = 10 , to minimize t)ie effect 
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Remark 5 

Note that many effective algorithms exist to find the eigenvalues of 
a matrix ;e.g. the QR family of algorithms); because of the large value of 
p occurring in M(P). it is essential however that balancing techniques [7] 
be used in finding the eigenvalues of the matrix M(p) to prevent problems 
of ill-conditioning arising in the eigenvalue calculation. Most existing 
eigenvalues packages have this feature (e.g. the EISPACK set of subroutines) 


5 . NUMERICAL EXAMPLE 

It is desired to find the transmission zeros of (C,A,B,D) where; 
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in this case, it is observed that r = m = 2 and that for a choice 
of X* = 2.17 (arbitrarily cUpI-.J [n-C..A-XM)-‘Bl - 2. so that the 

system is nondegenerate. On applying the algorithm of the previous 

section, the following eitoiivalues of (A + BK(-^ - DK) ^C) were obtained, 

for the case K ■■ P “1. on an Iu.l 370 digital computer using double precision 

arithmetic; tha eigenvalue al^oiithm used was a QR algorithm with balancing 
features (taken from the EISl’ACK set of subroutines). 
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Eigenvalues of A BK (^ - DK)~^C1 


p 

Eigenvalues 

»-> 

o 

00 

0.9999999899 

0.3411639019 ± jl. 161541399 
-0.6823278038 

-0.4999999949 ± jO. 1000000004 x 10^ 


0.9999999999 

0.3411639019 ± jl. 161541399 

• 

-0.6823278038 

-0.4999999999 ± jO. 1000000000 10^ 

loi'’ 

* 0.9999999999 

* 0.3411639019 ± jl. 161541399 

* -0.6823278038 

-0.4999999999 ± jO. 1000000000 ^ lo® 


*finite eigenvalues of A + - DK) 


In this case, it can be seen that: the transmission zeros of (C,A,B,D) 
are given by (to 10 significant figures): 1.000000000, 

0.3411639019 ± jl. 161541399, -0.6823278038. The transmission zeros obtained 
in this case agree with the exact values (given by the solution of 
(X-1) (X‘^+X+1) = 0) to at least 12 significant figures. 

6. CONCLUSION'S 

This paper has shown that for a general linear multivariable system 
described by (1), the poles of the closed loop system approach the trans- 
mission zero? of the system hhen high gain output feedback is used (theorem 2). 
Tliis result has formed the basis of a new algorithm for tlie numerical 
calculation of the transmission zeros of a general system. In the proposed 
algoritlim, the typical computation time/storage space required corresponds 
to that of finding the eigenvalues of a single n x n matrix, where n is 
the order of the system. The algorithm has been successfully used on systems 
of UJ) tu bill.. ul’OLl’. 
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APPENDIX 1 
Simplification of Equation (8) 


It is desired to show that (8) can be simplified to (9). It is as.«^umed _ 
that n=r and that (I^ - DKo) . K are nonsingular matrices. Let A. £ R 

b. c r"\ i = l,2,...,m, ® be defined as follows: 


A. ^ 
1 


' 0 1 0 ... 0 ' 


0 

0 0 1 ... 0 

A 

0 

• * * * 

• • • • 

• • • * 


• 

0 0 0 ... 1 


0 

0 0 0 ... 0 , 


. ^ . 


, .ID ^ Kp(I^ - DKo)’^C (la) 


- 1 . 


Then (8) can be written as: 

dot [A BKp(Ij. - DKp)"*C - Al] = 0 

which can be written as: 


dot block diag(Aj - Al /\^, - Al) 

•► block diag(b^ b^^l .H j = “ 


(2a) 


(3a) 
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which can be simplified to: 
[ n 

det 


-diag(X 

* 0 block diag 


m. 

♦ 

' 1 1 


1 1 

- 


X 


X 



• 


• 



n.'-l 

$ • • • 9 

n'-l 


_ 

X 1 

k J 


[x" J 



= 0 


(4a) 


which can be written as: 
det 


fl I D 111 

-diag(X “) + (p“^K~ -D)"^C 


block diag 




"in 


n'-l 
. m 


* 0 


(5a) 


Equation (5a) now can be written as: 


det 


-l.-l 


1 


( p'*K'*-D)‘-diag(X \...,X ®) + 


C block diag 



1 ' 


1 1 1 



X 


X 



• 


• 



• 

9 ••• f 

• 





n'-l 

X 

L'' J 



= 0 


(6a) 


or alternately as: 
det 


-p K diag(X "*) ♦ 


C block diag 




• •• •• 

♦ D diag(X "*) 


which can be written as: 
det 


n*-l 

. m 




= 0 


n. n 

-diag(X \...,X "*) ♦ pK 


n n 

♦ D diag(X ,...,X "’) 


C block diag 

« 0 


(7a) 


** 

r > 

1 


f y 

1 

• 


X 

• 


X 


_ 


» • • • » 

n'-l 

X 

W J 



or as 


n n 

dot (p K r(X) - diag(X ^....X "')] = 0 


(8a) 

(9a) 


which is the desired result. 
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ABSTRACT 

This paper considers the problem of observing a set of linear functionals of 
states for systems with unmeasurable disturbances. A transfer function approach 
is adopted to characterize the set of observers for a given system. The result 
also has direct application to state estimation in decentralized control systems, 

1. INTRODUCTION 

The problem of observing the state of a linear time-invariant multivariable 
system has been investigated by Luenberger in his early papers [1], [2]. Luen- 
berger has considered the design of minimal-order observers which estimate the 
whole state as well as the minimal-order observers which only estimate a given 
linear functional of the state. The problem of constructing a minimal-order 
observer which observes a given set of linear functionals of the state has been 
investigated in [3] - [6]. In [7] - [9]. different approaches are used to construct 
observers which can estimate the states of systems with unmeasurable disturbances. 
This paper considers the general problem of observing a given set of linear 
functionals of the state for systems with unmeasurable disturbances. In section II. 
a transfer function approach is adopted to characterize observers for a given 
system see [10] [11]. In Section III, it is shown that based on this transfer function 
characterization, one can then apply the minimization algorithm in [11] - [15] to 
obtain a lower bound for the order of dynamic observers required. For relatively 
low order systems, one can usually construct such observers with arbitrary poles 



without difficulty. However, the general problem of constructing minimal-order 
observers with arbitrary poles for systems with unmeasurable disturbances has 
not been fully resolved. In [6] Tarski's decision method has been applied to the 
design of minimal-order observers. For more details, see [16] , [17]. The 
result in this paper also has direct application to state estimation in decentralized 
control systems [18]. 

II. A TRANSFER FUNCTION CHARACTERIZATION OF OBSERVERS 
In this section we adopt the transfer function approach to the design of 
dynamic observers. This approach has been used by Retallack [10]. In the 
following. Theorem 1 gives a complete characterization of dynamic observers 
which observes a specified set of linear functionals of the state in terms of 
transfer function matrices for systems with unmeasurable inputs. 

Theorem 1 

Consider a completely controllable linear system specified by 

x(t) = Ax(t) + Bu(t) (1) 

y(t) = Cx(t) (2) 

where u(t) c is the input, x(t) c P*^ is the state, y(t) c P^ is the output and A , 

B and C are real constant matrices of appropriate size. 

Let A = diag [x..) be a mxm diagonal matrix with X^. = I if ujt) is measurable 

and X . = 0 if u.(t) is unmeasurable. Define u(t) = Au(t). 

11 1 

Consider another linear system with inputs u(t) c P*^ and y(t) e P^. ®nd output 
w(t) c P*^. 


2 



z(t) = Fz(t) + Gu(t) + Jy(t) 


(3) 


w(t) = Lz(t) + Mu(t) + Ny(t) (4) 

where z(t) c r" is the state, and F, G, J, L, M and N are real nnatrices of 

appropriate size, and the pair (L,F) is completely observable. (See Figure 1). 

Let K be an arbitrary real matrix of size qxn. Then the 
following two conditions are equivalent* 

(I) There exists a real number y ^ 0 and a map g; f” x f" -♦ [0. •) with 
g(O.O) = 0, such that Vx(o) c R*^ z(o) c r’^ • sind for any piecewise 

continuous control u; [ 0 ,») -♦ F^, we have 

Ilw(t) - Kx(t)|[ s: g(x(0), z(0)) X exp (-yt) Vt a 0. 

( II) The transfer function from u to w is equal to 

H 2 (s)Hj(s) = K(sl - A)"^B . (5) 

where 

H^(s) = 

H 2 ( 8 ) = L(sl - F)"^ [J [G] + [N |M] 
and 

Re X(F)s! -y if X(F) is a simple root of the minimal polynomial 

of F 

Re X(F) < -y if X(F) is a multiple root of the minimal polynomial 

ol F. 


■ 1 . 


C(sl - A) B 
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Remark 


Theorem 1 states that the class of linear multivariable systems with proper 
rational transfer functions, whose outputs w(t) asymptotically approach Kx(t) 
as t-* • , is characterized by the conditions of (II). 

In order to prove Theorem 1 , we first state the following lemma. 

Lemma 1 

Consider a completely observable system specified by 

x(t) = Ax(t) 
y(t) = Cx(t) 

where x(t) c R" is the state, y(t) c is the output, and A, C are real constant 
matrices of appropriate size. The following two conditions are equivalent; 

(i) There exists a function g; r” -* [0,*) with g(0) = 0, and a real number 
Y > 0 , such that 

I|C exp (At). x^|] s g (x^) exp ( -yt), Vx^cP". Vt a 0. 

(ii) Re X(A) s -y if X(A) is a simple root of the minimal polynomi.ial 

of A, 

Re X(A) < -y if X(A) is a multiple root of the minimal polynomial 

of A. 

The simple proof of Lemma 1 is omitted. 

Proof of Theorem 1 

(I) - (II) 

Let x(0) = 0 and z(0) r 0, then for any input u(. ), we have w(t) - Kx(t). 
Hence the transfer function from u to w is equal to K(sl - A) B. On setting 
x(0) s 0 and u(t) s 0 Yt i 0. 
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then 


w(t) = L exp (Ft) z(0) 

and 

x(t) = 0 Vt 2 0, 

Hence 

|Iw(t)II = ||L exp (Ft)x(0)[] s g (0. z(0)) exp (-Yt).Vt s 0. 

n * 

Since (L,F) is completely observable and z(0) e R is arbitrary, thenfrom 
Lemma 1 the result follows. 

(ID « { I) 

Let x^O) = 0 and z(0) = 0. Then since the transfer function from u to w 
is K(sl - A) we have 

0 = w(t) - Kx(t) Vt 2 0. (6) 

From (4), (6), and u = Au, 

0 = Lz(t) + MAu.t) + (NC - K) x(t) . (7) 

Assume that u( . ) has a discontinuity at time t and that both x( • ) and z( . ) are 
continuous functions at time t. In order to satisfy (7) , MA must be a zero 

matrix. 

Taking derivatives of (7). 

0 = LF z(t) + [LGA + NCB - KBj u(t) + [LJC + NCA - KA]x(t) (8) 

Using the same arguments as above, we conclude that the matrix [LGA+NCB- KB] 
0. Continuing this process, we have 

0 = LF^ ^z(t) + W.x(t), (1=1,2 n') (9) 

i 
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where W = NC-K, W = LJC + NCA-KA, etc. 

X M 

Eq. (9) can be written as 

Qz(t) = Wx(t) , (10) 



From the observability property of (L,F),the matrix Q Q is nonsingular. 

Hence from (10), there follows 

z(t) = (Q^Q)'^ Q^Wx(t) 

= Tx(t), Yt a 0, (11) 

where (Q^Q)’^Q^W. 

Next we calculate 

^ [z(t) - Tx(t)] = F[z(t) - Tx(t)] + [JC - TA + FT]x(t) + [GA - TB]u(t) (12) 
Assume that x(0) = 0 and z(0) = 0, from (11) and (12) , 

[JC - TA + FTjx(t) + [GA - TB]u(t) =0 Vt a 0. 

If u(. ) has a discontinuity at t, x(. ) is a continuous function at time t, hence 
[GA - TB] = 0. Furthermore, from the controllability property of (A,B), 
x(t) c p” can be any vector at time t > 0. Therefore [JC - TA + FT ] = 0. F rom 
(12), there follows 
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( 13 ) 


-^[z(t) - Tx(t)] = F[z(t) - Tx(t)] 
for any x(0), z(0), u(. ) and t ^ 0. 

Consider the following, 
w(t) - Kx(t) = Lz(t) + NCx(t) - Kx(t) 

= L[z(t) - Tx(t)] + [LT + NC - Kjx(t) 

= L-exp(Ft) [z(0) - Tx(0)] + [LT + NC - K]x(t). (14) 

From zero initial conditions x(0) = 0 and z(0) = 0, we know that 
w(t) - Kx(t) =0 Vt a 0. Then from the controllability property of (A,B), it 
follows that LT + NC - K = 0. Hence from (14), there exists a function 
g: r" x -» [ 0 , 08 ) with g(0,0) = 0 such that 

|lw(t)-Kx(t)lI = I{L exp(Ft) [z(0 )-Tx(0 )][j s g (x(0), z(0)) exp (-yt), 
n n' 

for all x(0) € P , z(0) c P and u(. ), and t a 0, where y is as specified in 
condition (II) of Theorem 1. Q. E. D. 

Remark ; Part of the above proof was motivated by results in [4]. 

III. LOWER BOUND ON THE ORDER OF OBSERVERS 
The minimization algorithm in [11] - [15] solves the following problem: let 

P(s) and Q(s) be given proper rational matrices of appropriate size: how does 
one find a proper rational matrix H(s), (assuming it exists), which satisfies the 
following equation 

P(s)H(s) = Q(s) 

such that H(s) is of least possible order (McMillan degree)? 

In the following, we will apply the minimization algorithm to find a lower 
bound on the order of a dynamic observer which is required to estimate a 
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specified set of linear functionals of the state of systems with unmeasurable 
inputs. 

From (5), we have 

C(sl - A)'^B 

A 

where C(sl - A) is the transfer function of the given system in (1) and (2), 
H 2 (s) is the transfer function of a dynamic system whose output tends to Kx(t), 
It is clear that if we take transpose of (5'), 




T T -1 T • T 

[B (si- A ) C I AjH^is) 


we can apply the minimization algorithm to the above equation to construct the 
T 

set of solutions H^Cs) of least McMillan degree. Moreover, as pointed out in 
[ 11 ], [ 12 ] the set of observability indices of all minimal-order solutions 1 ^ 2 ( 3 ) 
are equal. Since the solutions H^is) constructed by the minimization algorithm 
may not be stable only a lower bound on the order of the dynamic observer 
required is in general obtained. In the special case, however, when q = 1 or 
p = 1, one can rederive from (15) previous results on minimal-order observers 


[1] - [4]. For the general case q> 1 and p> 1, the problem of finding minimal- 
order dynamic observers with arbitrary pole locations is not yet fully resolved. 
As indicated in [6], Tarski's decision method seems to be suitable for this class 


of problems. For more details, see [16], [17]. 


IV. CONCLUSION 



required is derived. The result in this paper has direct applicatiPn to state 
estimation in decentralized control systems [18]. ^ 
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ABSTRACT 

This paper considers the problem of stabilizing a linear time-invariant 
multivariable system by using local feedback controllers and some limited 
information exchange among local stations. The problem of achieving a given 
degree of stability with minimum transmission cost is solved. 

1. Introduction 

When control theory is applied to solve problems of large systems, 
e.g., electric power systems, socioeconomic systems, etc., an important 
feature called decentralization often arises. Such systems have several 
local control stations; at each station, the controller observes only local 
system outputs and controls only local inputs. All the controllers are 
involved, however, in controlling the same large system. It is obviously 
important to know under what conditions there exists a set of appropriate 
local feedback control laws that will stabilize the complete system. Several 
different versions of this problem have been formulated and examined in 
[1], A necessary and sufficient condition for the existence of local output 
feedback control laws with dynamic compensation to stabilize a given system 
is given in [2]. This condition is stated in terms of a new notion, called 
"fixed modes" of a decentralized control system. If some element of the 
set of fixed modes is unstable, then the system cannot be stabilized via 
local controllers only. Some transmission of output information among local 

This work has been supported by the National Aeronautics and Space 
Administration, Langley Research Center under Grant No. NSG. 12i3 and the 
National Research Council of Canada, under Grant No. A4396. 
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stations is necessary. An algorithm wlucli minimizes the amount of trans- 


mission cost is given in Section 4. Section 

2 defines the problem 

and 

Section 3 summarizes some previous results. 

An example is given in 


Section .S. 



2. Problem Statement 



Consider a linear time-invariant multivariable system with v 

local 

control stations described by 



u 

x(t) = Ax(t) ♦ 

B.u. (t) 
1 1 ^ 

(la) 

rt 

II 

n 

X 

rr 

II 

l,...,o), 

(lb) 


1 . n ni ; n . 

whore x(t)e k is the state, u. (t)e R ^ and y.(t)c U ^ are the input 

and output, respectively, of the ith local control station (i = l,...,v). 

The matrices A, I^, and C^, li = l,...,v), are real, constant, and of 

appropriate size. I ho problem is to find v local output feedback control 

laws with dynamic cnirnensation in order to stabilize the whole system. The 

set of local feedback laws are assumed to be generated by the following 

feedback control lers: 


z.(t) = S.z.(t) Rjy.(t) (2a) 

u.(t) = Q.z.(t) ♦ K.y^(t) ♦ v.(t), 

<.,i = l,...,v,) (2b) 

n • m ^ 

wh<;rc z^Ult R ‘ is the state of the itii feedback controller, v^(t)c R ^ 

is tlie ith local external input, and .S^, R^, Q. , and K. are real constant 

matrices of appropriate size. Ii<|uatiun (2) can be written more compactly 


as follows: 
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ift) = Sz(t) ♦ Uy(t) (3a) 

u(t) = Qz(t) ♦ Ky(t) ♦ v(t). (3b) 

where 

S = block ilian 
U - block diiiK , . . . , 11^1 
Q ^ block cling iQj , • . . ,Q^1 
K ^ block diag (K ] 

z^t) ^ (Zj'‘'(t) ••• zj{t)] 

y^ft) = lyj‘^(t) ••• y^^(t)] 

u^(t) - [Uj^(t) ••• U^^(t)] 

and 

v'^(t) 4 (Vj'^u) ••• vj(t)). (4) 

When the feedback control law generated by (3) is applied to system (1), 
the closed-loop system is described by 


I 

X(t) 


A ♦ BKC BqI 

x(t) 



_z(t) 

1 

1 

z(t) 

+ j 

! 

! V(t). 

L J 


where 



V 


1 


B = 



m 


1 


m 

V 


and 



(5) 


(b) 


The set of local feedback control laws of (3) is t*^ be chosen so that 
the overall system (5) is asymptotically stable, (i.c., all poles of system 
(5) are in the open left-half complex plane Q ). This then is tlic problem of 
stabilizing a decentralized control system via local mil put feedback with 
dynamic compensation. 
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3. Solution to t h e Stabilization Problem - 

The following doflnitions are required to state the necessary and 
sufficient condition for the existence of solutions to the stabilization 
problem. 

Definition 1 

Consider the triple (C,A,B) e Rp*" X r"*" X r"^'" and the two sets of 
integers and 

m = L. and p = Z. ,p. 

i=l i * 1=1*1 

which specifies system (1). Let K be the set of block diagonal matrices 
as follows: 

= {K|K = block diag 
m xp. 

K^eR , i = l,**-,v). (7) 

Then the greatest cojnmon divisor of the sot of characteristic polynomials 
of (A 'IKC) , for all K , is called the fixed polynomial of the triple 
(C,\,B) with respect to.k", denoted by i;j(X;C,A, B,. fQ , i.e., 

•KA;C,A,B,K) = gcd{det(Xl - A - BKC)}. 

KcK 

Comment: llio map ii'( \;C,A.B,-) :jC*^gcd{dot (Xl - A - BKC)} is well defined 

nixD 

for any subset^, c R " and not just for the sot K of block diagonal natrices 
in (7J* 

Def init ion _2 

For a given triple (i:,A,U) £ rP*” X X r'''“" and the given set iC of 


s 


block diagonal matrices in (7), the set of fixed modus of (C,A,B) with respect 
to is defined as follows: 

= Da(A ♦ BKCl 
KeK 

where o[A ♦ HKC| denotes the set of eigenvalues of (A + HKC). liqui valently , 
A(C,A,S,X) can l>e defined us follows: 

A(C,A,H,)^) = {A|XeC and i|/(X;C,A,B,M = 0} . 

Theorem 1 [2| 

Consider the given system (C,A,B) of (1), with B arid C defined in (6). 
LetXbe the set of block diagonal matrices defined in (7). Then a nec- 
essary and sufficient condition for the existence of a set of local fee<lback 
control laws (3) such that the closed-loop system (5) is asymptotically 
stable is that 

A(C,A,B,K) C C 

where Q denotes the open left-half complex plane. 

Corollary 1 [2] 

Under the same assumptions as in Theorem 1, the necessary and suffic- 
ient condition fot the existence of a set of local feedback control laws 
in (3), such that all poles of the closed-loop system of (5) are in S 
is that 

A(C,A.B,fC)C S 

where S is any nonempty symmetric open subset on the complex plane (^(i.e., 
if X e S then its complex conjugate X* c S). 

For the proof of theorem 1 and some further details, please refer to 

12 ]. 

4. Minimum Transmission Cost 

In the previous sections, the problem of stabilizing u linear system 
via local output feedback is investigated. In case tliat a given system 



cannot be stabilized via local output feedback then it is necessary to transmit 
some observations of local system outputs to other conciol stations in 
order to stabilize the whole system. It is natural to assiimc that there 


is a certain amount of transmission cost, say > 0 , of transmittinij 
a scalar time function from station i to station j per unit time. 

The problem is to desi}>n a control scheme which stabilize a given i/steiii 
with minimum transmission cost. 

Consider the following sot of matrices 


j((p...i.j = l,---,V) = 


k|k = block{K.j} 


K„.‘ 


LKv,!--:kvJ 


, K,. e R 


m. xp. 
1 ‘j 


rank K. . = 


j ‘Pa 


) 


f8) 


The set of matrices^* defined above corresponds to the set of constant 
output feedback laws u ^ Ky. Hach feedback law u = Ky consists of a set 
of (sub) feedback laws u. = which requires the transmission of a 

time function with components from station j to station i. The 
total transmission cost required for the iniplcmentaticu of 
Kej((p.^,i,j = l,--.,v) is 



Hence the problem is to find a set of integers p.^, U < o^. < min (m^.p^) 
in order to minimize 

C “’ij-*,-*" 

subject to A(C,A,B.K(P| .v))C C • 


I-or a given set of p^j (i,j»l, . . . ,v) , the condition of (9) can be 
checked as follows: Uach submatrix K., of rank p. . can be written as 


7 


the product L..M.. where L. . is of size m.xp.. and M. . is 
' ij ij* ij I »j U 

Pj,jXp^ . Let f'C the set of cigciivaliics of A 

real parts. Then condition (i>) is satisfied if and only if 


n dct(>..I - A - B(block{L. .M. .) ICj 
ic{u 6} 


of size 

with iiOnnc-gative 


(to) 


is not identically zero. 


Remark 

The above condition implies that all unstable eigenvalues of A have been 
shifted by the applicatiorof feedback u = [ block } ] y, i. e. , none of unstable 
eigenvalues of A is a fixed mode. 

Since there is only a finite number of possibil ii i i-s in choosing L^j's, 
the minimization problem of (9) can be solved in a finite number of steps. 

Once the set (P^ . , i , j = l . . • • ,'») is determined, the actual construction 
of dynamic compensators is basically the same as in tlic proof ot il.corcm 1. 


5. Example 

Consider the following system with two control stations, 

f-Ql 


■’^1 ■ 


ol 

o 




“1 0" 

X.. 


0 0 0 


X, 


0 1 

2 




2 

+ 


X, 


0 0-2 




0 0 

3 




L 3 . 


L. -J 


u 


u 


1 

2 

1 . 


' 1 
“2 


y { =[0 0 1 ] 


”i " 


1 

0 

0 

2 

• 




^2 


0 

1 

0 _ 



3 


1 

uf 


and yj are the input and output at station 1; u^ and 


where 

and output ut station 2. Assume that the transmission cost or 


1 

yz 

2 


are the input 


11 = “22 = “• 


“l2 “ ' “21 ' 

In order to minimize the transmission cost, first consider the decent- 
ralized feedback 


- 1 , ^21 ^22 ~ 


= 


K 


k^l 0 0 


kzi 0 


' ^l'^21*^32’^3" ^ ^ 


^32 ^33 


With the above set K , the corresponding fixed polynomial of (C, A, B) is 

= g.c.d. 

W'32-’'33"' 

= > . 



X 

0 

-kii 


< det 

0 

;*^32 

X 

“^33 

1 

r 

, 1 

J 


This shows that with local feedback, there is a fixed mode at the origin. 
Hence the system can not be stabilized by local feedback. Now consider the 
case in which we transmit the output at station 1 to stn*’'nn 2. The corresponding 
transmission cost is 


^^^11 ^21 ■ ^22 °12 
And the set of feedback matrices is 

’f'“ir»2r“22='- “l2=0> 
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^1 

0 

0 


K 

K = 

*^21 

*^31 

0 

^2 

0 

. kj^, k^j, k^j, k^^. H > 


It can be shown that with respect to the above set of X. * the fixed 
polynomial of (C,A,B) is \. Hence, again the system can not be stabilized 
even with transmission of output from station I to station 2. 

Next consider transmitting a single linear combination of the output 

« 

at station 2 to station 1. The transmission cost is 


= ‘I’ll = '’u = ‘>21 = »21=“I = 2 


The corresponding set of feedback matrices is 



= Pl 2 = 

^21 = 

1, 0 

21 =°' 





'• 











^1 

^12 

^13 




s 

= i 

K 

K = 

^21 

^22 

^^23 

, all k..'s G B and rank 
ij 

^12 

i. 

^13 

k 

/ 

= 1 i 




0 

*^32 

^^33_ 

1 

^22 


‘ > 


The additional rank constraint on the matrix 
expressed as 


■'u ■'13 

■'22 *'23 


can be 


•■u ''13 

''22 ''23 


II 

Bl 


' TD 
IS) 

L. 


[v. v;i 


where Y, and Y are arbitrary real numbers, 

X w X w 

The corresponding fixed polynominal is 


g. c. d. 

k. ,'s 
ij 


f 

i det 

X -SiY, 
-®2''l 

■^1^2 ’*^11 
^■^2 ■‘^21 


-''32 

X■^2 


\ = 1 


REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 
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Hence there is no fixed mode in the system. By choosing 




= 0, and Y 2 ”^ 32 “ resulting system characteristic 


3 2 

polynomial is \ + 3\ + 3\ +1. Hence all three system poles have been 


assigned at -1. 


6. Cone 1 us ions 

In this pajici', the prol)lcm of stabilizing a large scale decentralized 
system via local controllers and with minimum transmission cost for information 
exchange among local stations is discussed. Clearly, a more efficient 
algorithm is required if one would like to apply this method to practical 
design problems. Several related publications are included in the reference. 
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